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ABSTRACT 

The structure of barotropically unstable disturbances in 
the tropics was examined with a two- level quasi-geostrophic 
model. The forecast equations were linearized and non- 
dimensionalized and the most unstable mode was found numeri- 
cally by use of the initial value technique. The wind profile 
Uj = a jV ec h 2 Cy/L) vvas use< ^* Simple CISK type heating was 
introduced to determine its effects on the growth of the 
waves. The jet wind profile was more unstable with easterly 
flow than westerly flow when no heating or friction was pres- 
ent. When the jet profile was present only at the upper 
level, the ratio of the disturbance amplitude at the lower 
level to the disturbance amplitude at the upper level decreased 
toward the equator. When the heating rate was not large enough 
to increase the growth rate, the amplitude ratio remained small. 
As the heating rate becomes larger, the ratio increases to a 
value larger than 1. 
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TABLE OF SYMBOLS 5 ABBREVIATIONS 



A e 

Bo 

C P 

f o 

g 



Vertical eddy viscosity 

Derivative of corcolis parameter at y 

Specific heat at constant pressure 

Coriolis parameter at y = 0 

Gravity 

dp/dt 



0 



V gz/ f Q 

R Gas constant for dry air 
p Density 

T s Average temperature from the standard atmosphere 
a (R 2 T s /p 2 g) (9T s /9z + g/C p ) 

C V 2 V 



z Height 

Q Heating added per unit mass 
V 4 Bi-harmonic operator (V 2 ) 2 
a Wind inflow angle 
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I. INTRODUCTION 



In the last few years there has been growing evidence 
that barotropic instability may be a common occurrence at 
the 200-mb level in the tropics. Krishnamurti [1970, 1971] 
observed a peak at wave numbers 6-8 in the wind spectrum at 
this level which he proposed was caused by barotropic insta- 
bility. A subsequent study by Kanamitsu, et al . [1972] showed 

that for the region between 15S and 15N waves of this scale 
do receive energy through barotropic interaction with the 
zonal mean flow and zonal wave number 1. 

Kuo [1949] examined the barotropic stability character- 
istics of a parallel zonal current and found that instability 
is possible only if the gradient of the absolute vorticity 
changes sign somewhere in the region. The 200-mb level in 
the tropics is the outflow level of the giant monsoon circu- 
lation which is induced by differential heating (both sensible 
and latent) due to land-sea distribution. The basic flow 
pattern thus varies in the planetary scale. Detailed obser- 
vations have found that synoptic- scale moving disturbances 
at this level tend to form in the easterly jet south of the 
Tibetan high [Krishnamurti, 1970] during the northern hemis- 
phere summer. This is a region of large vorticity gradient 
where the necessary condition for instability is occasionally 
satisfied locally. If these waves were the result of baro- 
tropic instability, they would extract energy from the mean 



8 



flow and the long waves, since the long waves and the mean 
flow combine to give large vorticity gradients in these 
regions. The energy exchanges computed by Kanamitsu et al. 
[ 1972 ] are consistent with this view. 

In this thesis the barotropic instability of an upper 
level easterly jet was examined. In addition to the calcu- 
lation of the growth rates, the downward propagation of energy 
was also investigated to examine the possibility that lower 
troposphere tropical disturbances may originate at the upper 
levels. Simple CISK type heating was introduced to ascertain 
how it affects the growth of the waves. 

This study uses a quasi-geostrophic two-level model. The 
equations are linearized and solved by the initial value 
technique. Growth rate and wave structure are extracted after 
an exponential growth is established. 
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II . THE FORECAST EQUATIONS 



The quasi-geostrophic prediction equations were used as 
the forecast equations. It is recognized that the quasi- 
geostrophic approximation is not very accurate near the 
Equator, but Matsuno [1966] has shown with a barotropic model 
that quasi-geostrophic motions are identifiable near the 
Equator if the wavelength is not too long. Also, in a study 
with a baroclinic model, Yamasaki [1969] has found that the 
error is small at twenty degrees latitude. 

The quasi-geostrophic equations are simplified through 
the use of a two-level model. The basic model is very similar 
to the models used by Williams et al. [1971] and Robertson 
[1973] . The equations were linearized and solved numerically 
by the initial value technique. The growth rates and wave 
structures were then computed. 

The simple two- level model was constructed by dividing 
the atmosphere into four layers with a constant pressure dif- 
ferential of Ap/2 (Fig. 1), numbered 0 to 4 from top to bottom. 
The boundary condition at p = p was to = 0 , which prohibits 
vertical propagation of energy. Charney [1969] has shown that 
such energy propagation is unlikely. At p = p Q the boundary 
condition was to = to where to was the frictionally induced to 
at the top of the Ekman layer. Charney and Eliassen [1949] 
derived the following expression for to g . 
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where A e is the kinematic eddy viscosity, a is the surface 
inflow angle, and is the surface geostrophic vorticity. 
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Fig. 1. Two-level prediction model 



It has been shown by Holton et al . [1971] and Chang [1973] 

that this expression is reasonably valid as far equatorward 
as ten degrees, but may be invalid close to the Equator south 
of the ITCZ. The surface geostrophic vorticity was approximated 
by 

C4 = Ca = V 2 T 3 (2.2) 

The quas i- geos trophic vorticity equation in pressure coordin- 
ates is 



|*j. V 2 T + Ikxvt • V(V 2 V) + 3 0 -|f - £ o -|| = o (2.3) 

where f Q = 2fisin ip Q , = 2 ft cos ^ Q /a, to = dp/dt, T = stream 
function, ip o is the reference latitude, and a is the radius 
of the earth. Observe that when to = 0 , (2.3) becomes the 
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barotropic vorticity equation which is sufficient to describe 
the barotropic instability. 

Applying this equation at levels 1 and 3 gives 



V 2 ^ + tKXVYj • V (V 2 'i 1 ) + e Q -|^ 






_ w O 

" f o ip = °’ 



(2.4) 



V 2 T, + (KXVT, • V(V 2 <iJ + 3 - f " M2 ^ = 0. 

dt 3 r 3 3 ' O dX 0 A P (2 5) 

If it is desired to stop at some intermediate point, say the 
tropopause, rather than to include the entire atmosphere, 
smaller layers could be used. 

Next, consider the quasi-geostrophic first law of thermo- 
dynamics in the form 



J)_ 3T 
3 1 3p 



[K X V T • V(^) + -js— to = 0 

dp I'o 



(2.6) 



where 






(2.6a) 



and T s is the temperature obtained from the standard atmosphere. 

The heating follows the formulation of Ooyama [1964] , 

Charney and Eliassen [1964], and Kuo [1965], i.e., the con- 
densation heating is proportional to the convergence of moisture 
in the Ekman layer. The form used was given by Ogura [1964]: 



Qa 



!. Cp T s 3 6 
2 e" 3p n w 4 



(2.7) 



where n is a non-dimensional parameter. This relationship was 
used in both rising and subsiding air representing a Fourier expansion 
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of the heating function. The factor 1/2 is therefore needed 
to maintain the correct amplitude of the heating perturbation. 
When (2.6) is appled at level 2 the result is 



" ^s) + IK X V 1 2 ^ * V ^i - - AP - f^ -- 2 - = 0. (2.8) 



Define the following quantities 



= 



+ V s 



y = - y 3 



which implies that 



= T M + T T ’ 






(2.9) 

( 2 . 10 ) 



(2.11) 

( 2 . 12 ) 



Here is proportional to the layer thickness and is therefore 
a measure of the mean temperature in the layer. Using these 
definitions, add (2.4) and (2.5) and divide the result by two, 
obtaining 

v2t m + I KXv ' } 'm ‘ v Cv 2 + |kxvt t ♦ v(v 2 t t ) 



+ 



ax 



f pu>e _ 
2Ap 



(2.13) 



For the second forecast equation, subtract (2.5) from (2.4) 
and eliminate w 2 by using (2.8) obtaining 
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- u 2 ) T t 



+ |K X VV M * V(V 2 - y 2 )^ + |K X VT - ) 




(2.14) 



where 



2 




(2.14a) 



These are the prediction equations for the model. 

Although the normal mode method gives the phase speeds, growth 
rates, and wave structures for all modes, the initial value 
approach gives the phase speed, growth rate, and wave structure 
of only the most unstable mode. This approach was chosen 
because only information about the most unstable mode was 
needed. 

It is convenient to write ^ and in the following form: 

= ^ (^ + My,t) cos kx + B(y,t) sin kx , (2.15) 

= F(y) + C(y,t) cos kx + D(y,t) sin kx. (2.16) 

where k is the x wave number. The coefficients A through D 
are the Fourier amplitudes of the disturbance and E and F are 
the mean zonal fields. 

These expressions for and ipj are substituted into 
(2.13). The various sine and cosine terms are then separated, 
and all products of the quantities A through D neglected. 
Equating the coefficients of the cosine terms gives 
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JL filj - Ak 2 ! 
3 1 ''TTy r AJCJ 



3 E 3 2 P . 3 F 3 2 D 



= k[^ 0 -; 



3 y 3 y^<- 3 y 3 y" 



a 3 E p 3 j_F 

3 y 3 D 3 y 3 u 



,2 f 3 E R + 3 F 

‘ k ( sy B + sy 



D )] - B o Bk - K %* 



- k 2 ) (A - C). 
( 2 . 17 ) 



For the sine terms the result is 



3 f 3 2 B 
3 1 l 3y2 



- Bk 2 )= k [- 



3 E 3 2 A 3 F 3 2 C + l!l A , 3 3 F 



3 y 3 y 2 3 y 3 y 2 3 y 3 



3 y : 



+ k 2 ( 9E 



k-A + a — C)] + 3 Ak - K(^ 5 - 
3 y 3 y J J o v 3 y 2 



- k 2 ) (B - D) 
( 2 . 18 ) 



Repeat the procedure for equation ( 2 . 14 ) for the cosine terms 
which gives 



Jt (0 - ck2 - cv2 l m k tli fly + f£ |1| - t§D U 2 + P 2 ) 



3 y 3 y 2 3 y 3 y 2 3 y 



3JF 

3y 



B (k 2 - U 2 ) 



3 3 E 
3 y 3 



D 



3 3 F 

3 y 3 



B ] 



- 8 0 Dk + (K - S) (|y 2 - k 2 ) (A - C) . 



( 2 . 19 ) 



The results for the sine terms are 



3 , 3 2 D 



(’Ey - Dk2 - dp 2 ) = k [ffc(k 2 + p 2 ) + -^a (k 2 - p 2 ) 



3 F 



3 1 3 y 



3 y 



3 y 



3 E 3 2 C 3 F 3 2 A . 3 3 F „ , 3 3 E 

+ - — - A + . _ ■ , L J 



3 y 3 y 2 



3 y 3 y 2 3 y( 



3 y • 



+ 3 0 Ck + (K - S)(|-- - k 2 ) (B - D), 



where 



K = £qS.. AV2 

2RTc 1 f J 



( 2 . 20 ) 



( 2 . 20 a) 
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and 



S = ■§■ 



30 Po 



A m V2 



0 



3p 



(“T~) 
0 o 



( 2 . 21 ) 



The computational equation for vertical velocity is given 
below : 



00 9 



3 ljjrr. 

[3T- 1 + 



3iJ/ 



T 



Apo L 3t M 3x 



+ v Sl|, T + KQ^ 
v m 3y 2£ 0 J > 



( 2 , 22 ) 



where 



u 



m 



v 



m 



1 

f 0 



1 

f 0 3x 



(2.23) 

(2. 24) 
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Ill . BOUNDARY CONDITIONS AND FINITE 
DIFFERENCING SCHEME 



The finite difference scheme is illustrated below with 
a sample variable M: 



9M 

ay 




M 



i- 1 



) , 



(3.1) 



9 2 M 
9y2 




2M . + M. ,) 

l l- 1 J 



(3.2) 



9 3 M 
9y 3 



1 

2H 3 



K M i + ; 



2M i + i 



+ M.) 

1 J 



(M. - 2M. , + M. ,,)] 

v i l-l i- 2' J 

(3.3) 



where i is the grid index and H is the distance between grid 
points . 

Centered time differences were used for all quantities 
except those involving friction and heating. The friction 
and heating terms were computed at time (t - At) . The first 
step in all cases is a forward time step. The second order 
equations for time tendencies were solved by the exact method 
of Richtmyer [1957, p. 101], 

Yanai and Nitta [1968a] studied the problem of using 
finite difference equations to solve dynamic instability 
problems of non-divergent barotropic currents with boundaries. 
In this study rigid boundaries are placed at y = ± W/2 where 
W is the total width of the computational field. The boundary 
conditions were A=B=C=D=Oaty=±W/2. 
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Fig. 2. Mean wind profile Uj = a jU Q sech 2 (y/L) 
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IV. WIND PROFILE AND INITIAL CONDITIONS 



This study uses the following zonal wind profile: 

U. = a.U sech 2 (y/L) , (4.1) 

■} j o v ' 

where a^ can be ± 1 or 0. The barotropic stability of this 
profile has been analyzed by Lipps [1962]. Figure 2 contains 
the profile (4.1) for a^ = 1. 

The initial conditions are given below: 



A = sin iry/W 




(4.2) 


B = 0 




(4.3) 


C = 0 




(4.4) 


D = 0 




(4.5) 


E = a i + a s £ u L tanh 

2 oo 


(y/L) 


(4.6) 


F = a i - a 

-i— * f U L tanh 

2 oo 


(y/L) 


(4.7) 



Note that the E and F fields are independent of time. 
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V. NON- DIMENSIONALI ZATION PROCEDURE 



In order to simplify the interpretation of the results, 
the equations are written in non-dimensional form as follows 
x = Lx^, y = Ly'', and t = L/t x . The stream function is 
written 

¥ m = " LU 0 /u M Cr) d r + ^jCx^y^t") (5.1) 

and 

= - LU o /u t (y '') dy ^ (5.2) 



when these expressions are substituted into (2.13) and 
(2.14), and when the products of the ip's are dropped, the 
non-dimensional equations can be written: 




u 4 )v2 *m 



+ 



(3 



d 2 u . 
dy 2 ' 




TV 2 (ip. 



M 




0 (5.3) 



^3t + U 3x^ v2 " e ^T + ^ " dy 2 -* 3x" + C<3 - r)V 2 ( ip M * iJ; T ) 

(5.4) 



where 



* - o ~ i n foL sin 2n A p 

e - LV, b = B 0 L 2 U . r - — 



(5.5) 



and 



6 = Umax L p 0 g d 2 1 gp I s ^ n ^ a A e ^ 2 

U4foT ^ 2f ^ 

u o o A 4 o 



(5.6) 



The quantity e is a rotational Froude number, 3 is a non- 
dimensional form of 3, T is the ratio of the square root of 
the Ekman number [Greenspan, 1968] to the Rossby number, and 
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6 is the heating parameter. The primes on the independent 
variables have been dropped. These equations are equivalent 
to the set 2.17, 2.18, 2.19, and 2.20 which were solved 
numerically . 
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VI. COMPUTATIONAL PROCEDURE 



The following constants were used in the computations: 



o 



o 



= 100 centibars, 

J 

L = 400 kilometers, 

W = 4000 kilometers, 

U_ = 10 meters per second, 

= varied with latitude, 

= 3.0 meters per second squared per 
centibar squared, 

= varied with latitude, 

= 50 centibars, 

= 10 meters per second squared, 

= 22.5 degrees, 

= varied with latitude 
= 50 kilometers, 

= 1.0 hours , 

= 4, 6, 8 
3¥/9p= -0.6 degrees per centibar. 

A forecast period of 100 days was used so that the most unstable 
mode would have reached the maximum growth rate. For conven- 
ience, growth rates were computed over the last two days of 
the forecast period. This was done by assuming the amplitudes 
could be written as 

Xj = a Q exp (n U q /L) t 1 , (6.1) 



3 

Ap 

A e 

a 

p 2 

H 

At 

n 



max 



x , = a 0 exp ( n U ^ /L ) t o> 



o 



( 6 . 2 ) 
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where n is 
By forming 



the non-dimensional growth rate and a Q is a constant, 
the ratio \ 2 /x l , the growth rate can be shown to be 

L In (xa/xj) (( 

n " u 0 '(tf- tj- ' (6 ' 3) 
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VII. RESULTS 



A series of experiments were conducted over a wide variety 
of conditions. The more important experiments are given in 
Table I as a function of aj , a 3 , § , c, 6, T, n, and k", where 
n represents the maximum growth rate for each set of conditions 
and k' is. the corresponding non-dimensional wave number kL. 

When some of these parameters were computed over a range of 
values, the range is given. 

Experiment 1 was conducted to determine the difference in 
growth rate characteristics between easterly and westerly jets. 
In this experiment no heating or friction was included and the 
mean flow was the same at both levels. In order to obtain a 
continuous variation between the easterly and westerly flows, 
the non-dimensional parameter g (see equation 5.5) was defined 
to assume the sign of U. Thus f3 is positive for westerly flow 
and negative for easterly flow. In the remaining experiments 
the jet will be easterly, but g will be regarded as positive. 
For each g the maximum growth rate and the corresponding wave 
number were determined. Figure 3 contains these quantities 
as a function of g; the maximum growth rate of n = 0.1888 and 
occurs at g = -0.352. The corresponding wave number is k" = 
0.99. This value of k^ is close to the minimum k” which occurs 
when 3 = -0.112. Fig. 3 shows clearly that the easterly jet 
is more unstable than the westerly jet. In particular, the 
growth rates for the easterlies over a range of beta are 
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greater than the value for beta = 0. This is in contrast to 
the normal notion that the presence of beta has a stabilizing 
effect for instabilities. 

Experiment 2 investigated the downward propagation of 
energy from a barotropically unstable disturbance as a func- 
tion of the latitudinal location of the jet. An easterly jet 
was introduced at the upper level and the lower level mean 
flow was set to zero. Surface friction was included, but 
heating was neglected. A measure of the downward propagation 
is the ratio R = |Y,| / I ¥, I where these amplitudes are 

obtained from the exponentially unstable solution. Figure 4 
contains isolines of R as a function of k" and latitude. The 
curve n shows the wave number which has the maximum growth 
rate at each latitude. The R value decreases as one moves 
from high to low latitude as predicted by Charney [1963] . 

The R value also decreases with increasing wave number at 
each latitude. This figure shows the dependence on epsilon. 

Experiments 3, 4, and 5 all contain heating and were con- 
ducted with a fixed latitude of 10 degrees, otherwise they 
are the same as experiment 2. In these experiments a weighting 
factor was applied to the heating following Bates [1970] and 
Yamasaki and IVada [1972], The following expression was used 

n = n max /cosh 2 y/L. (7.1) 

This form implies that moisture is available for the heating 
in the region | y | < L. If this was not done, CISK heating would 
occur in regions away from the basic jet. In experiment 3 
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where n =4 the non-dimensional heating parameter is 
‘max 

6 = 0.165. Figure 5 shows the growth rate (n) and the ampli- 
tude ratio (R) as a function of the non-dimensional wave number 
k ' . Note that the maximum growth rate of n = 0.187 occurs at 
k' = 1.0. These two values are very close to the values ob- 
tained without heating. The R is fairly large at the lower 
values of k"' with a maximum at 0.5. R becomes very small in 
the area of maximum growth. Apparently the barotropic insta- 
bility effects dominate in the areas of maximum growth while 
the heating is relatively more important for k'' £ 0.5 as 
evidenced by the larger values of R. This effect is even 
more pronounced in experiment 4. 

In experiment 4 where n =6 the non-dimensional heating 
^ ‘max ° 

parameter is $ = 0.247. Figure 6 shows that the maximum growth 
rate of n = 0.186 occurs at k^ = 1.0 just as in experiment 3. 

As k^ varies away from 1.0 to either extreme the growth rate 
approaches 0.167. The R is larger in both of these regions 
where the largest value occurs for the largest value of k'. 

A secondary maximum in R occurs at k' = 0.6. This indicates 
that heating dominates the solution in the regions which have 
the larger values of R. 

In experiment 5 where p = 8 the non-dimensional heating 
parameter is 6 = 0.330. Figure 7 shows that the growth rate 
is fairly constant over the k' considered. This is similar 
to the results obtained by Williams and Robertson [1973] who 
acquired essentially a constant growth rate with a hyperbolic 
tangent wind profile using larger values of 6. 
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GROWTH RATE n 




Fig. 3. The growth rate n (solid line) and the corre- 
sponding wave number (dashed line) as functions 
of 3 for experiment 1. Westerly flow occurs for 
3 > 0 and easterly flow for 3 < 0. 
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DEGREES LATITUDE 



EPSILON 




Fig. 4. Isolines of R = lT,| /It.I (curved 
b 1 3 i max 1 i 'max v 

lines) as functions of k", latitude, and 

epsilon . 
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Fig. 5. Growth rate (dashed line) and amplitude ratio 
(solid line) as a function of wave number for 
experiment 3. 
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GROWTH RATE n 




Fig. 6. Growth rate (dashed line) and amplitude ratio 
(solid line) as a function of wave number for 
experiment 4. 
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Fig. 7. Growth rate (dashed line) and amplitude ratio 
(solid line) as a function of wave number for 
experiment 5. 
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VIII. CONCLUSIONS 



Recent studies suggest that disturbances may form through 
barotropic instability in the easterly jet south of the 
Tibetan High at 200-mb. In this thesis the instability of an 
analytical jet profile was investigated with a two- level 
quasi -geostrophic model. The downward propagation of energy 
from a barotropically unstable upper level disturbance was 
investigated to see if it could initiate lower tropospheric 
disturbances . 

The hyperbolic secant squared wind profile was used in 
all of the numerical experiments. The first experiment, which 
contained no heating or friction, compared the growth rates 
with easterly and westerly flow. The easterly jet was the 
more unstable and its maximum growth rate was well above the 
growth rate for g = 0. The energy propagation was studied 
by introducing the jet at the upper level with no mean flow 
below. The ratio (R) of the disturbance amplitude at the 
lower level to the disturbance amplitude at the upper level 
was computed. With surface friction included, R decreased 
from high latitudes to low latitudes; at each latitude R 
decreased with increasing wave number. When CISK type heating 
was introduced with a large enough heating coefficient, the 
growth rate became independent of wave number. In this case 
R had a finite value and it increased with increasing wave 
number. For smaller values of the heating coefficient the 
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barotropic instability controls the growth rate and R remains 
very small. 

The result that the easterly jet is more unstable than 
the westerly jet supports the view that barotropic instabil- 
ity occurs south of the Tibetan High. The small value of the 
amplitude ratio at low latitudes agrees with Charney's [1963] 
conclusion that the vertical energy propagation is small. 

The inclusion of CISK heating did increase the low level 
amplitude, but this could not occur unless a finite amplitude 
were present at the lower level. Since these results were 
obtained with a quasi-geostrophic model with only two levels, 
it is important to carry out further studies using more 
levels and also allowing for gravity waves. This would 
require a multi-level primitive equation model in either 
spherical coordinates or on an equatorial 6-plane. A more 
sophisticated heating parameterization, such as that by 
Arakawa and Schubert [1974], should be used. 
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